Ap femtoscopy in collisions of Ar+KCl at 1J6A GeV 
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Results on Ap femtoscopy are reported at the lowest energy so far. At a beam energy of 1.76 A GeV, the 
reaction Ar+KCl was studied with the High Acceptance Di-Electron Spectrometer (HADES) at SIS18/GSI. A 
high-statistics and high-purity A sample was collected, allowing for the investigation of Ap correlations at small 
relative momenta. The experimental correlation function is compared to corresponding model calculations al- 
lowing the determination of the space-time extent of the Ap emission source. The Ap source radius is found 
slightly smaller than the pp correlation radius for a similar collision system. The present Ap radius is signif- 
icantly smaller than that found for Au+Au/Pb+Pb collisions in the AGS, SPS and RHIC energy domains, but 
larger than that observed for electroproduction from He. Taking into account all available data, we find the Ap 
source radius to increase almost linearly with the number of participants to the power of one-third. 



PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Gz 

Two-particle correlations of hadrons at small relative mo- 
menta are widely used to study the space-time extent of their 
source created in heavy-ion collisions or other reactions in- 
volving hadrons (for review see [1 1). Generally, the sign and 
strength of the correlation is affected by i) the strong interac- 
tion, ii) the Coulomb interaction if charged particles are in- 
volved, and iii) the quantum statistics in case of identical par- 
ticles (Pauli suppression for fermions, Bose-Einstein enhance- 



ment for bosons). Their interplay makes, e.g., pp correlation 
functions rather complex with a suppression at very small rel- 
ative momenta (due to items ii) and iii)) followed by a max- 
imum which results from the short-range attractive potential 
of their strong interaction. In contrast, correlations of non- 
identical hadrons, with at least one partner being uncharged, 
are sensitive to the strong interaction alone. Hence, Ap cor- 
relation functions are well suited to study the spatiotemporal 
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extension of the particle emitting source, provided the Ap in- 
teraction is known j2). Presently, information on Ap corre- 
lations is rather scarce. Due to the necessity of having high 
statistics of protons and A's at small relative momenta, such 
information is available at high beam energies only. Experi- 
mental Ap correlation functions have been reported for central 
collisions of Au+Au/Pb+Pb by E895 at the Alternating 
Gradient Synchrotron (AGS) at Brookhaven National Labo- 
ratory (BNL), by NA49 ||4) at the Super Proton Synchrotron 
(SPS) at CERN, by STAR at the Relativistic Heavy Ion 
Collider (RHIC) at BNL, and, for electroproduction from He, 
by CLAS [6 1 at the Thomas Jefferson National Accelerator 
Facility (JLab). On the other hand, the higher the energy the 
wider the particles are distributed in momentum space result- 
ing in a reduced probability of finding pairs with small rel- 
ative momenta. Thus, at lower beam energies an exception- 
ally large number of events and the more compact momentum 
distribution of A hyperons might compensate partially their 
lower production probability. Indeed, such a high-statistics 
A sample is available to us Q ID making an analysis of Ap 
correlations feasible. Furthermore, at low energies the feed- 
down correction to the A is smaller and more reliable, being 
of importance for comparisons to models applying explicitly 
to primary A's like the model by Lednicky and Lyuboshitz [9 1 
which we use in the present investigation. 

In this paper we report on the first observation of Ap corre- 
lations in heavy-ion collisions at beam energies below 2 GeV 
per nucleon. The experiment was performed with the High 
Acceptance Di-Electron Spectrometer (HADES) [10] at the 
Schwerionensynchrotron SIS 18 at GSI, Darmstadt. HADES, 
primarily designed to measure di-electrons JJTJ, offers also 
excellent hadron identification capabilities 1171 l8l fT2l [131 . 

A 40 Ar beam of about 10 6 particles/s with kinetic energy 
of 1.756A GeV (^snn — 2.61 GeV) was incident on a four- 
fold segmented target of natural KC1 with a total thickness of 
5 mm corresponding to 3.3 % interaction probability. The data 
readout was started by a first-level trigger (LVL1) decision, 
requiring a minimum charged-particle multiplicity > 16 in the 
time-of-flight detectors covering polar angles from 18 to 85 
degrees. The integrated cross section selected by this trigger 
comprises approximately the most central 35 % of the total 
reaction cross section. The corresponding mean number of 
participants, estimated with the UrQMD transport approach 
Ifl4l amounts to (A part ) = 38.5 ± 2.5. About 700 million 
LVL1 events were processed for the present investigation. 

Let Fi2(Pi, P2) be the coincidence yield of pairs of parti- 
cles having the momenta pi and P2. Then the two-particle 
correlation function is defined as 

s-,/ \ . r ^-^events, pairs ^12 (Pi ) P2 ) 

C(pi,p 2 ) =A/ ^ — 7 r. (1) 

2^ events, pairs * 12, mix (Pi ) P 2 ) 

The sum runs over all events and over all pairs satisfying 
certain conditions. Event mixing, denoted by the subscript 
"mix", means to take particle 1 and particle 2 from different 
events belonging to the same multiplicity interval of 8 units 



width [12]. Due to limited statistics, the six-dimensional cor- 
relation function is projected onto one-half of the relative mo- 
mentum in the pair cm. frame, i.e. we consider C(k) with 
k = |pi — P2I/2. The normalization factor AT in ([TJ is 
fixed by the requirement C(fc) — > 1 at large relative momenta, 
k = 0.10 - 0.15GeV/c. The statistical errors of C(k) are gov- 
erned by those of the coincidence yield, since the mixed yield 
can be generated with much higher statistics. 
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FIG. 1: Top: The pn~ invariant mass distribution. Light and dark 
hatched areas represent the statistics of the A signal and the corre- 
sponding A impurity, respectively, both contributing to the Ap cor- 
relation analysis. Bottom: The uncorrected Ap correlation function. 



In the present analysis we identified the A hyperons through 
their decay A — > p7r~ . To allow for a clear A selection vari- 
ous cuts were applied (cf. [8]). The resulting invariant mass 
distribution of events containing at least one additional well- 
identified proton is displayed in the upper panel of Fig.[TJ 

Taking the resulting A sample (hatched area, ±6 MeV in- 
terval around A peak), we started the correlation study by 
combining, for each event containing a A candidate, the A 
with those protons not already contributing to the A. The re- 
sult is a Ap relative-momentum distribution containing a total 
of 240,000 proton-A pairs. However, only 2,700 (260) pairs 
contribute to the interesting region of small relative momenta, 
k < 0.1 (0.04) GeV/c. The resulting raw Ap correlation func- 
tion is displayed in the lower panel of Fig.[TJ A clear enhance- 
ment at small relative momenta is found, indicating the effect 
of the strong Ap interaction. 

Corrections for the finite small-angle acceptance losses 
were deduced from simulations. 25 million thermal A's (one 
per event) were generated using the fireball option of the event 
generator PLUTO 03]. The fireball parameters of the A are 
chosen such that the simulation reproduces both the experi- 
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FIG. 2: The same as Fig.[TJ but for simulated A hyperons with their 
decay products embedded into true experimental data. The full curve 
in the lower panel is a fit with a Fermi-type function (see text) used 
to correct for reconstruction losses (due to close tracks) of the exper- 
imental Ap correlation function displayed in Fig.[TJ 



mental values of the effective inverse slope parameter at mid- 
rapidity and the rapidity width ll7l[8l [T2ll . The simulation data 
are processed through GEANT 1161 . modeling the detector 
response. The GEANT data were embedded into real experi- 
mental data and processed through the full analysis chain. The 
correlation function of, ab initio, uncorrelated pairs of (simu- 
lated) A hyperons and (experimental) protons is expected to 
be flat at unity. Any deviation from this value has to be in- 
terpreted as detector bias due to finite granularity, momen- 
tum and near-track resolutions. Figure[2]shows the correlation 
function of simulated A's combined with protons from exper- 
imental data. A slight suppression at small relative momenta 
is visible indicating the bias of the apparatus and the track- 
finding losses. The full curve is a fit with a Fermi-type func- 
tion, facc(k) = 1 + Ai/(1 + e'^^)^ 3 ), used to correct the 
experimental Ap correlation function for reconstruction losses 
due to close tracks. The close-track correction is performed 
via C(k) — > C(k)/ f acc (k), while the purity correction is done 
by the transformation C(fc) — > 1 + (C(fc) — 1) / Pair Purity. 
Here, PairPurity is the product of the proton (0.95 ± 0.02) 
lfT2l and A purities. The A purity, in turn, follows from the 
product of the signal purity of 0.816 ± 0.013 as derived from 
the p7r~ invariant mass distribution of Fig.[T](which increases 
to 0.876±0.014 for half the given mass window) and the feed- 
down correction of the decay E° — > A7. (Feed down from 
decays of heavier particles is neglected in the present investi- 
gation.) On the one hand, the shape and the magnitude of the 
S°p and Ap correlation functions are predicted to be close to 
each other [17]. On the other hand, residual correlations car- 



ried by secondary A's are expected to be partially dissolved 
by decay kinematics. Hence, we approximated the feed-down 
correction midway between these boundaries, yielding a fac- 
tor of 0.90±0.10. Here, the upper limit of the correction factor 
(i.e. no correction) is related to identical Ap and E°p corre- 
lations combining into the measurable one, while the lower 
limit meets the case of maximum feed-down correction, i.e. 
assuming a complete loss of the S°p correlation after the elec- 
tromagnetic decay of the S°. The corresponding yield ratio 
of A/(A + E°) of about 0.80 was derived from calculations 
with the UrQMD transport approach lfl4l . At last, utilizing 
GEANT simulations, the influence of the momentum resolu- 
tion onto the Ap correlation function was investigated. It ap- 
pears to be negligible (i.e. less than 2 % increase of the width 
of C(k) — 1, which can be approximated by a Gaussian func- 
tion with a dispersion of (27 ± 3) MeV/c), as expected from 
the narrow widths of both the A and S~ signal peaks lf8l . 

Finally, we applied a fit to the corrected Ap correlation 
function (Fig. [3] left panel) with a fit function derived from 
[9|. Following the procedure described in (3), we used the 
Ap scattering lengths (/q = — 2.88fm, /g = — 1.66fm) and 
effective ranges (dp = 2.92 fm, dy = 3.78 fm) for the spin sin- 
glet and triplet states of the Ap system as given in [2|. Notice 
that by increasing the source radius both the height and the 
width of the correlation above unity decrease. We applied the 
new log-likelihood minimization procedure recently proposed 
in lfT8l the result of which, however, deviates hardly from that 
of a x 2 minimization. The optimum Gaussian radius provided 
by this fit amounts to 

r = (2.09 ± 0.16 ifcg ig;3 123?) (2) 

where the 1st error is the statistical error, while the 2nd, 3rd, 
and 4th ones represent the systematic errors due to the uncer- 
tainties of the close-track correction with embedded A's, due 
to the pair purity correction (with a systematically smaller ra- 
dius by 0. 12 fm for half the A mass window as given in Fig.[T]i, 
and due to a ±25 % variation of the scattering lengths entering 
the model (5), respectively. (Note that a fit to the uncorrected 
correlation function as displayed in the lower panel of Fig.[T] 
would give a Gaussian source radius of 3.1 fm.) 

We studied the source radius resulting from the fit when 
changing the whole set of Ap singlet and triplet scattering 
lengths and effective ranges describing the Ap potential. In 
ref. |[T9l a compilation of such parameters is given. For the 
sets by the chiral effective field theory EFT '06, the Jiilich '04 
and Nijmegen NSC97f models the fit yields a radius of (2.13± 

0. 15) fm, (2.15 ± 0.15) fm and (2.11 ± 0.15) fm, respectively. 
Though there are quite substantial differences in the parame- 
ter sets, the fitted source radii do not show serious differences 
and are very similar to the value Q derived above for the set 
from Q. Obviously, any sensitivity to the peculiarities of the 
Ap two-body interaction is dissolved after the integration over 
the corresponding wave functions and the source distribution, 

1. e. the procedure yielding the theoretical correlation func- 
tions. Hence, the correlation observed in the Ap momentum 
difference is primarily sensitive to the source size. 
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FIG. 3: Left: The Ap correlation function after close-track and pu- 
rity corrections. The full curve represents the best fit (see text) with 
the Analytical Model by Lednicky and Lyuboshitz |9|. Right: The 
Gaussian radius of the Ap emission source as a function of system 
size. The symbols indicate data taken with HADES at SIS, CLAS 1 6 1 
at JLab, E895 Q] at AGS, NA49 g) at SPS, and STAR [5 1 at RHIC, 
respectively. The dashed line is a linear regression to the data. 



The Ap source radius Q appears slightly smaller than the 
pp source radii derived by the FOPI collaboration ||20l com- 
paring two-proton small-angle correlations in central colli- 
sions of Ca+Ca to the predictions of the Koonin model [21 1. 

The present Ap source radius may be compared to the cor- 
responding radii derived in other experiments. In the right 
panel of Fig. [3] we show the Gaussian radius as a function of 

1 /3 

the number of participants to the power of one-third, A part , 
which is calculated from the centrality and the total size of 
the corresponding collision system using a geometrical model 
of penetrating sharp spheres. While for the data measured 
by NA49 at SPS (158.4 GeV Pb+Pb, preliminary), by 
STAR at RHIC (Au+Au at y^vw = 200 GeV), and by 
CLAS [6 1 at JLab (preliminary results from e+ 3 He( 4 He) at 
4.7 (4.46) GeV), the Gaussian radius r is determined using 
the same model as in the present analysis, the half-maximum 
radius R 1/2 derived by E895 at AGS (6 A GeV Au+Au) ap- 
plying an imaging procedure was transformed to a Gaussian 
radius via r$ = R1/2/ V21n2. Clearly, the Ap source radius 
increases with system size. Similarly to the systematics of 
two-pion HI [22] and two-kaon [23 1 source radii, we find an 

1 /3 

almost linear increase with A part . Finally, we point out that 
the precision, i.e. the relative error, of our Ap source radius 
Q can well compete with or is even better than the accura- 
cies reported at higher energies. Nevertheless, with presently 
available statistics it is not possible to use the data sample for 
getting deeper insight into the specifics of the Ap interaction. 

In summary, we observed, for the first time in heavy-ion 
collisions at SIS energies, Ap correlations at small relative 
momenta. The system Ar+KCl at 1.76 A GeV beam energy 
was measured with the HADES detector. The Ap correla- 
tion function, after proper corrections, was compared to the 
output of the model by Lednicky and Lyuboshitz. The result- 
ing Gaussian source radius is slightly smaller than that de- 
rived from pp correlations of a similar collision system, when 



interpreted with the Koonin model. The present Ap source 
radius is found significantly smaller than that deduced from 
Au+Au/Pb+Pb collisions at higher beam energies and larger 
than that reported for electroproduction from He. It increases 

1 /3 

almost with A part . A significant variation of the scattering 
lengths and the effective ranges of the Ap interaction entering 
the model changes only marginally the source radius result- 
ing from the one-parameter fit to the experimental Ap correla- 
tion function. This indicates the robustness of the correlation 
function w.r.t. uncertainties in the knowledge of the Ap two- 
body interaction. Finally, we note that the envisaged high- 
intensity runs of both the upgraded HADES and the Com- 
pressed Baryon Matter Experiment planned at the future Fa- 
cility for Antiproton and Ion Research (FAIR) [24 1 may pro- 
vide the statistics necessary for more differential correlation 
analyses even of pairs of rare particle species, hence gaining 
information not only on the extents of the emission sources 
but also on their shape or spatiotemporal separation. 
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